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1 Introduction 

The evolution of primordial black holes in scalar-tensor theories has been studied in the literature [Q, g, 
||, f| . The question of what happens to a black hole in an expanding universe in these theories was posed 
by Barrow, who discussed two scenarios (a) The effective gravitational "constant" G(t) at a black 
hole changes along with its cosmological evolution so that the size of a black hole is approximated by 
R = 2G(t)M . (b) G remains constant at the black hole event horizon while it evolves on larger scales; a 
large inhomogeneity in G is therefore generated. The case (b) was called "gravitational memory" because 
the black hole remembers the value of G at its formation time. In either case constraints on primordial 
black holes would be modified 

Later Jacobson claimed that there is no "gravitational memory" effect, analyzing the evolution of 
a scalar field 4>{t,r) in Schwarzschild background ||. His analytic solution showed that </> at the event 
horizon evolves along with its asymptotic value (f>(t, oo). However, the solution used has a special form 
with very fast variation of G oc t^ 1 in the background universe. It was also argued that even if the black 
hole mass in the Einstein frame is constant the mass in the Jordan frame is time-dependent. 

In order to investigate the time-dependence of the black hole mass, Saida and Soda constructed a 
"cell lattice" universe in the Brans-Dicke (BD) theory Q. In their model the universe are tessellated 
by identical polyhedrons, which are replaced by Schwarzschild-like black holes. It was shown that the 
black hole mass has adiabatic time dependence, and that its time dependence is qualitatively different 
according to the sign of the curvature of the universe. 

As an extension of Saida and Soda's work, we consider a "Swiss cheese" model j5j in the BD theory 
and discuss the evolution of the radius and mass of black holes. The ordinary Swiss cheese model refers 
to a cosmological model in which spherical regions in the Friedmann- Robertson- Walker (FRW) universe 
are replaced by Schwarzschild black spacetimes. Here we construct such a model in BD cosmology. 

2 Background Universe in Brans-Dicke Theory 

The BD theory is described by the action, 
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where <f> is the BD field, lj is the BD parameter, and C m is the matter Lagrangian. The variation of 
Eq.(|l]) with respect to and <fr yield the field equations: 



f 8?r uj 
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+ V^V^- <? M „D0, (2) 



D * - 2^3^ ® 
As a background universe we assume the FRW spacetime: 

ds 2 = -dt 2 + a(t) 2 l-^r + r 2 (d# 2 + sin 2 9d<p 2 )\ , (4) 
( 1 — kr z J 

where fc is the sign of the spatial curvature. As an energy-momentum tensor we consider a dust fluid: 

= pu^Uv, (5) 

where p and are the density and the four velocity of dust, respectively. Then the field equations (^|) 
and (|3|) reduce the following equations for the background universe: 

-*+?-5-*!+5($)' < 6 > 

where an overdot denotes (i/ei/j and H = a/a is the Hubble parameter. 

In the case of A; = we know the analytic solution for the dust universe: 

a(t) = a (t-t+) x +(t-t_) x -, 4>{t) = Mt-t+Y + {t-t-T-, (8) 

with 

uo + 1 ± Jl + 2w/3 l±3v/l + 2w/3 

A± = ' K± = s^Tl ' (9) 

where ao, <po and t± are arbitrary constants. If we take t+ = t_, the general solution (||) reduces to the 
special power-law solution, 

a(t) = a (t - t )%£* , 0(t) =0o(t-to)5&i. (10) 

In the limit of t — > oo, the general solution (|J) converges to the special power-law solution ([n]). If 
the present cosmic age t is large enough, observations cannot constrain the relation between t + and i_. 
Therefore, we keep t + — t_ a free parameter. 

In vacuum case (p — 0) we have the analytic solutions for all k. The vacuum solution for k = is 
expressed as 

a(t)=a t^J, = oi Tfe, (11) 

with 

1 ± VI + 2^/3 

a = > (12) 

where we have omitted the arbitrary constant to by fixing the origin of the time coordinate t. Introducing 
the conformal time rj — J dt/a, the vacuum solutions for k = ±1 are expressed as 

k = +1 : a (77) = (sin 77) (cos 77) , ^(77) = (tan 77) A , (13) 

k = — 1 : 0(77) = (sinh 77) (cosh 77) , ^(77) = (tanh?7) A , (14) 

with 

i = ± 3W (15) 
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3 Modified Swiss-Cheese Model 



We now consider the model of a black hole embedded in the FRW universe. We replace a sphere in the 
FRW spacetime with a vacuum region which contains a black hole. Here "vacuum" means T^ v — 0, and 
does not imply that IZ^ = due to the existence of the BD field. 

Extending Israel's junction conditions for a singular (or regular) hypersurface ||, Sakai and Maeda 
have studied bubble dynamics in the inflationary universe 0. It was found that one can solve the 
equations of motion for the boundary without knowing the interior metric if the interior is vacuum, 
T^v = 0, or has only vacuum energy (a cosmological constant), T M „ = —pg^u- Applying this method to 
the present model, we find the mass and the radius of a black hole without specifying the interior metric, 
as we shall show below. 

Let us consider a spherical hypersurface S which divides a spacetime into two regions, V + (outside) 
and V~ (inside). We define a unit space-like vector, N^, which is orthogonal to £ and points from V~ 
to V + . To describe the behavior of the boundary, we introduce a Gaussian normal coordinate system, 
(n, x l ) = (n, r, 9, <p), where r is chosen to be the proper time on the boundary. Hereafter, we denote by 
^E' ± the value of any field variable $ defined on £ by taking limits from V 4 . 

For the matter field, we consider dust (or vacuum) for V + and vacuum for V~: 



T + 



o. 



(16) 



Although we assume a smooth boundary at which there is no surface density, it is not obvious that this 
matching is possible at all times. Therefore, we introduce the surface energy-momentum tensor on the 
boundary surface, 



lira 



dn Tij — diag(— a, w, w) 



(17) 



where a and w are the surface energy density and the surface pressure of £, respectively. 

Introducing the extrinsic curvature tensor of the world hypersurface £, Kij 
junction conditions on £ as H 



Ni-j, we can write the 



[Kij] 41 — — I Sij 



2uj + 3 



(18) 



2oj + 3 



(19) 
(20) 



where we have defined the bracket as [^ r ] ± = ^ + — W and the vertical bar | as the three-dimensional 
covariant derivative. The junction condition for the BD field is derived from Eq.(||) as 



24tt 
"3 + 2lo 



TrS*, 



(21) 



which implies that <f> is continuous at £ and inhomogeneous in 

The extrinsic curvature tensor of £ in the homogeneous side V + is expressed as |7| 



e = <y(l + vHR) 
9 R 




— R? \ H 2 



(22) 
(23) 



where 



R = a(t)r\s, v 



dr 
i— 
dt 



a7 



and e = sign(i^g) = sign 



OR 
dn 



(24) 
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From Eqs.@, @, (|T|), @, @-( ||), we obtain the equations of motion: 

dR dr „ . 

— = —v + HR, (25) 
at ax 

f /_ „ \ „ 2w*l 27TCT f (l-2w) 2 ] 7V0 



dt '{\ J R dx J <t> I (2w + 3) 

do- 2cr(l + w) dR . . 

where w = vo/a. 

Once initial values of i?, u, and a are given, the equations of motion (p5|)-(^7|) determine their evolution. 
As discussed in Ref. , initial values should satisfy the angular component of (|l^) , 



n , „ m _ L{dR\ 2 R MS STraR m + l + w 

7(1 + ui/i?) - e f + UJ ~^r = - 2c+3 ' (28) 

where we choose e~ = +1. Rms is defined as 

RMS^R-il-g^R-R-), (29) 

where i?~ is defined as R~ = ^Jgee at £ on the side of V - . Because the Misner- Sharp mass is defined as 
I 

M MS =^(l-g^ = R-=R-), (30) 

we call Rms the "Misner-Sharp radius" . Note that Rms is a purely geometrical quantity and independent 
of theories of gravitation. 

If we considered a bubble in which there is no black hole (or singularity), we would have to solve 
the field equations with regularity condition at the centre and the boundary condition (pg)), as done 
in Ref.0. However, because we are interested in black hole solutions, we do not have to take such a 
regularity condition into account. Thus we can use Eq.(|2^) to determine Rms- 

At the initial time we suppose v = and a = zu = 0. Then Rms is expressed as 

Rms = R 3 (h 2 + ^ . (31) 

Let us now discuss whether v and er remain zero during the ensuing evolution. Suppose w — 0, then 
the only nontrivial term in Eq.(|26|) is r y 2 v 2 p/<r. If v and a evolved from zero, Eq.(p8|) shows vH ~ (r/4>, 
so that pv 2 /(T r~j pv/Htp. Therefore, Eqs.(|2q) and fl27| ) guarantee that, if v = a = at a certain time, 
v = a = at all time. Interestingly, this result is true only for the dust case, w/a = 0; otherwise 
(2w/R)(dr/dx) in Eq.(p6|) would shift v from zero. 

In the case of Schwarzschild spacetime, the Misner-Sharp radius coincides with the event horizon. 
Although it is not always true for general spacetimes, we speculate that the Misner-Sharp radius is a 
well-defined measure of the size of a black hole. In the next section, we show the evolution of Rms f° r 
black holes in several background cosmological models. 

4 Evolution of Black Holes 

The evolution of the Misner-Sharp radius for the k = dust universe is given by Eqs.(§) and (|3ll), 

Rms = ay Q + 2 (t - t+) 3A + {t _ t _ )3 A_ (32) 

where ro is the comoving radius of the vacuum region. Equation ( |32"| ) shows that the black hole size 
decreases with time. 
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If we define the black hole mass as 

M MS = ^, (33) 
it coincides with the mass defined by Saida and Soda Q. For the k = dust universe, we obtain 

M MS = ^ + j^-)\t t + f^+(t *_)»>-+"-. (34) 

It is easy to see that Eq.(|34|) reduces to Mms — constant, if we choose t + = £_, which is the same result 
as that found by Saida and Soda Q. They also showed Mms increases for k = + 1 and decreases for 
k = — 1, and concluded that the evolution of the mass depends qualitatively on the sign of the curvature 
of the universe. We should note, however, that their conclusion is true only for the special case t+ = i_, 
or equivalently, only for the asymptotic behavior of Mms at t — > oo. 

Next, let us consider the scalar-field dominated (vacuum) universe. The evolution of Rms and Mms 
in flat, open, and closed universes are given by 

Q 

f 3-l-A \ | A 

k = +l: Rms = y( cosr ?) "(sin^) ~(cos2r? - sin2r? - A), (37) 

"^_|_^ 1-|-A 

Mms = -j{t&nri) x (cosr))~^-(smr)) ~ (cos It] - sin 2?7 - A) , (38) 

3 

fe = -l: i? A /s = ^■(coshr;) _1 ^(sinh7;)~ i ^ A (cosh277 -sinh27;- A), (39) 

7^ 3-l-A 1 -j- A 

Mms = ^■(tanh77) A (cosh?7) _ T _ (sinh?7) _ ~(cosh2ry - sinh2?7 - A). (40) 

If we take u> > 500, both a and A are negligible. We see that Rms and Mms both decrease with increasing 
time, except in the contracting phase of the k = +1 universe. 



5 Discussion 

We have constructed a modified "Swiss cheese" model in the Brans-Dicke theory, and discussed the 
evolution of black holes for dust and vacuum universes. We have defined the size of a black hole by the 
Misner-Sharp mass, and found that it always decreases as long as the universe is in an expanding phase. 

Although we have not specified the metric around a black hole, the mass and radius of a black hole 
we have obtained coincide with those of Saida and Soda Q, who assumed a Schwarzschild-like metric. 
This means that their ansatz of the Schwarzschild-like metric does not introduce a specialization of the 
problem. One fundamental problem remains. Although the radius Rms represents a typical size of a 
black hole, the relation between it and the event horizon has not been demonstrated explicitly and this 
problem will be considered further in future work. 
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